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Abstract 

Background: During September 2009, a large dust storm was experienced in Sydney, New South Wales, Australia. 
Extremely high levels of particulate matter were recorded, with daily average levels of coarse matter (<10 |jm) 
peaking over 1 1,000 Lig/m 3 and fine (<2.5 Lim) over 1,600 Lig/m 3 . We conducted an analysis to determine whether 
the dust storm was associated with increases in all-cause, cardiovascular, respiratory and asthma-related emergency 
department presentations and hospital admissions. 

Methods: We used distributed-lag Poisson generalized models to analyse the emergency department presentations 
and hospital admissions adjusted for pollutants, humidity, temperature and day of week and seasonal effects to 
obtain estimates of relative risks associated with the dust storm. 

Results: The dust storm period was associated with large increases in asthma emergency department visits (relative 
risk 1.23, 95% confidence interval 1.10-1.38, p < 0.01), and to a lesser extent, all emergency department visits 
(relative risk 1.04, 95% confidence interval 1.03-1.06, p<0.01) and respiratory emergency department visits (relative 
risk 1.20, 95% confidence interval 1.15-1.26, p < 0.01). There was no significant increase in cardiovascular emergency 
department visits (p = 0.09) or hospital admissions for any reason. Age-specific analyses showed the dust storm was 
associated with increases in all-cause and respiratory emergency department visits in the >65 year age group; the 
<5 year group had higher risks of all-cause, respiratory and asthma-related emergency department presentations. 

Conclusions: We recommend public health measures, especially targeting asthmatics, should be implemented 
during future dust storm events. 
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Background 

A large growing body of literature supports the associ- 
ation between exposure to particulate air pollution and 
adverse health outcomes [1-10]. Reports discriminate, to 
some extent, between health outcomes associated with 
different sources of particulate matter (PM). Relatively 
few published reports concentrate on exposure to primar- 
ily coarse respirable particles (PM 10 ), less than 10 um in 
diameter [4,10-14], These reports are generally associated 
with geological dust, in areas where dust exposure is fre- 
quently recurring and there are a large number of days of 
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moderate to high exposure available to investigate poten- 
tial health outcomes. 

An issue with many studies is the potential mitigating 
effect of health warnings on high pollution days. Such 
warnings are designed to affect population behaviour to 
reduce the incidence of adverse health outcomes. Warn- 
ings are commonly disseminated during high PM air 
pollution episodes. On Wednesday 23rd September 2009, 
the city of Sydney (in New South Wales (NSW), Australia), 
experienced a rare large dust storm [15]. Extremely high 
levels of PM were recorded, with levels of PM 10 and fine 
PM (less than 2.5 um in diameter, PM 2 .s) peaking at over 
11,000 |-ig/m 3 and 1,600 (ig/m 3 (daily average levels), re- 
spectively. PM 10 and PM 2 .5 levels experienced during this 
dust storm exceeded any extreme PM levels during bush 
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fires and dust storms in the previous 15 years, by several 
orders of magnitude [7,8]. The dust storm dissipated con- 
siderably by early evening but two days later another less 
severe dust storm event lasting approximately 24 hours 
occurred. Exhaustive media coverage advised of health 
risks and mitigation options - avoid exposure if pos- 
sible; stay indoors (in air conditioned buildings); do not 
exercise; follow asthma plans and seek medical help if 
respiratory or cardiovascular symptoms occurred. These 
messages were widely disseminated in mainstream 
media, and SMS text alerts were sent to those subscrib- 
ing to the NSW Office of Environment and Heritage 
(OEH) air quality alert website [16]. 

Because of the sudden onset of extreme levels of PM, 
this dust storm provides an opportunity to assess mor- 
bidities during an extreme air pollution event that re- 
ceived widespread media attention. We conducted an 
analysis to determine associations between this dust 
storm and outcomes of emergency department (ED) pre- 
sentations and hospital admissions for all causes and 
specific sub-groups of cardiovascular disease, respiratory 
disease and asthma. 

Methods 

Hospital and air pollution data 

Data were restricted to the Sydney Statistical Division 
(SSD) from 1 January 2004 to 31 December 2009 (a total 
of 2191 days). The SSD is a standard geographical clas- 
sification defined by the Australia Bureau of Statistics 
and contains the metropolitan area of Sydney (the lar- 
gest city in Australia and the capital of NSW). In 2008, 
the SSD had an estimated population of 4,339,722 
persons, consisting of 1,423,521 households with the 
majority of residents living in an urban environment 
(97.5%). 18.9% of the population were aged 14 and 
under and 12.1% aged 65 and over [17]. 

Daily data on ED presentations and hospitalisations 
were obtained from the NSW Ministry of Health (NSW 
Health) and included information from 29 public hospitals 
in the SSD (over 95% of ED presentations in this region, 
pers. comm. NSW Health). Planned or pre-arranged visits 
were excluded from the analyses. The data was classified 
into cardiovascular disease (ICD-9-CM: 036.42, 074.22, 
093.2, 098.84, 112.81, 115.04, 115.14, 115.94, 390-459, 
785, 786.5, 793.2, 794.3; ICD-10-AM: B37.6, G45, 100-199, 
R00-R03, R07.1-R07.4, R93.1, R94.3, Z86.7), respiratory 
disease (ICD-9-CM: 033, 460-519, 768-770, 786.0-786.4, 
799.1; ICD-10-AM: A37, J00-J99, P20-P28, R04.2-R06, 
R09.0-R09.3, R09.89) and asthma (ICD-9-CM: 493; ICD- 
10- AM: J45-J46) daily counts. 

Daily air pollution data was obtained from NSW OEH 
on PM 10 and PM 2 .5, carbon monoxide (CO), nitrogen di- 
oxide (N0 2 ), sulfur dioxide (S0 2 ), ozone (0 3 ), humidity 
and temperature from 19 air pollution monitoring stations 



located in the SSD. The pollutant data were summarised 
as daily means of data from monitoring stations during 
the study period. The pollutant and meteorological data 
from the recording stations provided by OEH were 
complete for the SSD area. 

In this study we used routinely collected de-identified 
administrative hospitalisation and ED visit datasets and 
therefore we did not require approval from an institu- 
tional ethics committee. The NSW Ministry of Health, 
as the data custodian for these two datasets, approved 
the use of the datasets. 

Statistical analysis 

Data was analysed using R 2.15.1 (The R Foundation for 
Statistical Computing) with the 'dlnin package. 

Descriptive statistics for air pollutants, humidity, 
temperature, ED presentations and hospital admissions 
for each of the 4 month periods during 1 August 2004 
to 30 November 2009 were analysed to compare with 
corresponding levels on the day of the dust storm. 
August is the final month of winter and November is 
the end of spring in the southern hemisphere, these are 
the surrounding months of the dust storm. Outcomes 
were total number of daily ED presentations and daily 
hospital admissions for all-causes and for cardiovascular 
disease, respiratory disease and asthma. 

An indicator variable from 23 September 2009 to 6 
October 2009 (the dust storm period) was created to es- 
timate associations between the dust storm and the ED 
presentations and hospital admissions. The high PM 
level on the day of the storm is problematic to model, 
the indicator variable is designed to capture the extreme 
PM levels occurring on the day of the storm and after- 
effects. CO, N0 2 , S0 2) 0 3 , humidity and temperature 
levels were considered as potential confounders. These 
air pollution and meteorological variables were included 
as distributed lags, where multiple lags are simultan- 
eously included in the time series model [18-20]. Dis- 
tributed lag models have been used to previously model 
effect of air pollutants on health effects [21-24]. An indi- 
cator variable for the day of week was included to cap- 
ture daily variation. To model seasonality and remove 
autocorrelation, we used a natural cubic spline smooth- 
ing function [25] for the date of the count. Degrees of 
freedom for the spline functions were selected via as- 
sessment of model fit using based on Akaike information 
criterion and residual autocorrelation. 

The statistical analyses used distributed-lag Poisson 
Generalized Models [26] to model the count data. The 
maximum distributed-lag and degrees of freedom for the 
cubic spline were chosen by assessing model fit with 
modified Akaike and Bayesian information criteria [23]. 
Relative risks (RR), 95% confidence intervals (CIs) and 
p-values (Wald Type III p-values) were calculated. The 
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results presented in this article are from the final model 
derived from this variable selection process, listings of 
variables included in the model are reported in tables in 
Additional file 1. 

We examined residual plots (quantile-quantile plots, 
histograms and plots of residuals against fitted and addi- 
tive predictors) and goodness of fit statistics to check 
model diagnostics [26]. We also examined the effect of 
the dust storm period on two age groups (0-5 years, 65 
years and older) and on sex. We performed sensitivity 
analyses on the model results by different smoothing 
methods and altering window length of the indicator for 
the dust storm period. 

Results 

Cardiovascular disease ED presentations comprised 8.4%, 
respiratory 9.8% and asthma 1.2% of total ED presenta- 
tions; cardiovascular disease admissions 9.9%, respiratory 
12.9% and asthma 1.2% of total hospital admissions. 
Descriptive statistics for daily mean pollutant levels, 
temperature and humidity are presented in Table 1, 
showing comparisons of concentrations on the dust 
storm day (23/9/09) with levels in the dust storm period 
(23/9/09-6/10/09), 2009 and previous years (2004- 
2009). PM levels on the day of the dust storm were ex- 
tremely high (daily average PM 10 = 11 705 (ig/m 3 ; daily 
average PM 2 . 5 = 1638 ug/m 3 ). A second day of elevated 
PM levels occurred on 26/9/09 (daily average PM 10 = 
783 (ig/m 3 ; daily average PM 2 .5 = 110 ug/m 3 ). The levels 
of PM during the dust storm were of an unprecedented 
order of magnitude higher than those experienced dur- 
ing previous years. Other pollutants, humidity and 
temperature levels were similar to levels in 2009 and 
previous years. There were no bushfire events, heat 
waves or other unusual metrological phenomena during 
the dust storm period [27]. Correlations between air 
pollutants, temperature and humidity are reported in 
Additional file 1: Table SI. 

The RRs for ED presentations from the modeling are 
shown in Table 2. The dust storm period is associated 
with a 4.3% increase in risk of all-cause ED presenta- 
tions, compared to the remainder of the study period. 
The dust storm period was not significantly associated 
with cardiovascular ED presentations, but was associated 
with a 20% risk increase in respiratory ED presentations 
and a 23% increased risk increase of asthma ED presen- 
tations compared to other periods. 

The corresponding RRs for hospital admissions from 
the modeling are shown in Table 2. The dust storm 
period was not associated with any significant increased 
risk of all-cause, cardiovascular or asthma hospital ad- 
missions. The dust storm period was associated with a 
10% decrease in risk of respiratory hospital admissions, 
compared to other periods. 



Results from the age-specific analyses for ED presenta- 
tions are presented below. Amongst the >65 years group, 
the RR for all-cause ED presentations were increased 
compared to the results for all ages (RR 1.109, 95% CI 
1.075-1.144, p < .001). There was no association between 
the dust storm and cardiovascular ED presentations (RR 
0.940, 95% CI 0.872-1.013, p = 0.107). There were signifi- 
cant associations between the dust storm period and 
respiratory ED presentations (RR 1.281, 95% CI 1.174- 
1.397, p < .001) and none associated with asthma ED 
presentations in the >65 year group (RR 1.362, 95% CI 
0.981-1.891, p = 0.065). 

In the <5 year age group the RRs for all-cause ED pre- 
sentations were increased compared to results from all 
ages and >65 year age group (RR 1.237, 95% CI 1.190- 
1.286, p < .001). RRs were not calculated for cardiovascular 
disease as there were too few events to calculate meaning- 
ful statistics. Compared to the >65 year age group, the <5 
year age group had similar risk of respiratory (RR 1.286, 
95% CI 1.192-1.392, p < .001) and asthma ED presenta- 
tions (RR 1.268, 95% CI 1.041-1.545, p = 0.019). 

Results from the sex-specific analyses for ED presenta- 
tions are presented below. Increased risks were associ- 
ated between the dust storm period and all-cause ED 
presentations for both sexes (compared to other pe- 
riods), these effects are similar between the sexes (males: 
RR 1.095, 95% CI 1.074-1.116, p < .001; females: RR 
1.144, 95% CI 1.121-1.168, p < .001). There were no sig- 
nificant associations between the dust storm period and 
cardiovascular ED presentations in either sex (males: RR 
0.976, 95% CI 0.912-1.045, p = 0.489; females: RR 0.988, 
95% CI 0.919-1.062, p = 0.749). The positive associations 
between the dust storm period and respiratory ED pre- 
sentations are also similar between the sexes (males: RR 
1.282, 95% CI 1.205-1.365, p < .001; females: RR 1.318, 
95% CI 1.233-1.409, p < .001). Asthma hospital admis- 
sions for both sexes showed significant increased risks 
during the dust storm period (males: RR 1.254, 95% CI 
1.066-1.475, p < .001; females: RR 1.418, 95% CI 1.227- 
1.638, p < .001), compared to other periods. 

The results remained consistent regardless of smooth- 
ing methods used (splines, loess curves) for seasonal ef- 
fects and distributed-lags, and to varying the length of 
the time window for the period of the dust storm. 

Discussion 

We conducted an analysis to investigate whether a se- 
vere dust storm period was associated with increases in 
ED presentations and hospital admissions. After control- 
ling for potential confounders, the dust storm period 
was associated with large significant increases in asthma 
ED presentations, and to a lesser extent, respiratory ED 
presentations, all-cause ED presentations. There were no 
significant increases in ED visits for cardiovascular 
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Table 1 Descriptive statistics of pollutant levels 



Pollutant 


N 


Min 


Median 


IQR 


Mean 3 


Std Dev 


Max 


August-November 2004-2008 


CO (ppm) 


610 


0.06 


0.36 


0.29 


0.42 


0.23 


1.6 


N0 2 (ppbhm) 


610 


0.55 


2.05 


0.79 


2.05 


0.61 


4.15 


0 3 1 hourly (ppbhm) 


610 


1.57 


3.38 


1.08 


3.62 


0.95 


8.73 


PM 10 (Mg/m 3 ) 


610 


10.57 


34.69 


18.53 


38.84 


23.13 


258.18 


PM Z5 (M9/m 3 ) 


610 


3.3 


13.83 


9.47 


16.08 


12.39 


1 74.42 


S0 2 (ppbhm) 


610 


0 


0.21 


0.21 


0.25 


0.17 


0.9 


Rel. humidity (%) 


610 


26.5 


65.89 


16.67 


64.63 


12.58 


91.26 


Temperature (C) 


610 


7.12 


15.45 


5.98 


15.66 


3.95 


28.28 


August-November 2009 


CO (ppm) 


122 


0.13 


0.47 


0.36 


0.54 


0.25 


1.24 


N0 2 (ppbhm) 


122 


0.25 


1.81 


0.92 


1.85 


0.62 


3.33 


0 3 1 hourly (ppbhm) 


122 


2.03 


3.25 


0.73 


3.42 


0.73 


5.99 


PM 10 (Mg/m 3 ) 


122 


13.95 


38.48 


31.42 


165.75 


1156.62 


1 1 704.8 


PM Z5 (Mg/m 3 ) 


122 


4.22 


14.66 


15.53 


34.13 


161.05 


1637.82 


S0 2 (ppbhm) 


122 


0 


0.22 


0.18 


0.24 


0.18 


0.89 


Rel. humidity (%) 


122 


35.29 


65.97 


18.4 


63.88 


14.16 


89 


Temperature (C) 


122 


8.37 


15.37 


5.05 


15.41 


3.41 


27.5 


23/9/09-6/10/09 


CO (ppm) 


14 


0.13 


0.36 


0.3 


0.44 


0.23 


0.82 


N0 2 (ppbhm) 


14 


0.25 


1.49 


1.06 


1.63 


0.83 


3.18 


0 3 1 hourly (ppbhm) 


14 


2.7 


3.11 


0.34 


3.2 


0.41 


4.12 


PM 10 (Mg/m 3 ) 


14 


14.07 


47.57 


43.75 


925.25 


3108.95 


1 1 704.8 


PM Z5 (Mg/m 3 ) 


14 


5.58 


12.5 


18.3 


137.85 


432.57 


1637.82 


S0 2 (ppbhm) 


14 


0 


0.11 


0.19 


0.13 


0.1 


0.29 


Rel. humidity (%) 


14 


35.58 


52.2 


38.12 


59.51 


20.64 


88.54 


Temperature (C) 


14 


12.2 


14.65 


3.27 


15.19 


2.46 


20.8 



23/09/2009 b 



CO (ppm) 1 


0.33 


N0 2 (ppbhm) 1 


0.56 


0 3 1 hourly (ppbhm) 1 


2.98 


PM 10 (Mg/m 3 ) 1 


1 1 704.8 


PM Z5 (Mg/m 3 ) 1 


1637.82 


S0 2 (ppbhm) 1 


0.06 


Rel. humidity (%) 1 


48 


Temperature (C) 1 


17.52 



Arithmetic mean. b Consists of 1 measurement for the day 23/09/09. 



disease. The dust storm period was not associated with 
an increase in asthma ED visits in the >65 year age 
group. It is possible that older people affected by asthma 
were more likely to have stayed indoors and provided 
with extra care. They were therefore less likely to be ad- 
versely affected by the dust storm. The dust storm was 
not associated with increases in risk for all-cause, cardio- 
vascular or asthma hospital admissions. It is possible 



issues were not severe enough for subsequent admission 
to hospitals. 

This dust storm originated in the Western NSW and 
South Australian arid zones and was clearly tracked over 
time during its travel path. There were no bush fire 
events or other unusual metrological phenomena during 
this period [27]. The study event is one very extreme 
event and it is difficult to find similar events to compare 
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Table 2 Adjusted relative risks of ED presentations and hospital admissions associated with the 2009 dust storm 



ED Presentations 



Hospital Admissions 



Category 


RR 


Lower 95% CI 
All-cause 


Upper 95% CI 


p-value 


Category 


RR 


Lower 95% CI 
All-cause 


Upper 95% CI 


p-value 


23/9/09-6/10/09 


1.043 


1.028 


1.058 


<.001 


23/9/09-6/10/09 


0.987 


0.970 


1.005 


0.149 


Non dust storm 


1.000 


Cardiovascular 






Non dust storm 


1.000 


Cardiovascular 






23/9/09-6/10/09 


0.957 


0.910 


1.007 


0.092 


23/9/09-6/10/09 


0.985 


0.936 


1.038 


0.581 


Non dust storm 


1.000 


Respiratory 






Non dust storm 


1.000 


Respiratory 






23/9/09-6/10/09 


1.199 


1.145 


1.255 


<.001 


23/9/09-6/10/09 


0.898 


0.853 


0.945 


<.001 


Non dust storm 


1.000 


Asthma 






Non dust storm 


1.000 


Asthma 






23/9/09-6/10/09 


1.230 


1.099 


1.377 


<.001 


23/9/09-6/10/09 


1.141 


0.991 


1.313 


0.066 


Non dust storm 


1.000 








Non dust storm 


1.000 









Details of variables included and distributed-lags are included in Additional file 1: Tables S2-S9. 



to. The Sydney region is prone to bushfires and these 
are best for comparison, although the levels of PM dur- 
ing these fires are 10 times lower. This dust storm 
appears to be the first well recorded event worldwide 
within a metropolitan area and therefore allows unique 
opportunities to study the health effects of the storm on 
the occupants of the metropolitan area [8]. 

Our results for all-cause, cardiovascular, and asthma 
hospital admissions are consistent with the current 
literature. Other studies have failed to find significant as- 
sociations between dust storms and cardiovascular ED 
presentations and hospital admissions [28,29], which 
agree with our findings. We found increases in respira- 
tory and asthma ED presentations during the dust storm 
period, which is in agreement with similar studies 
[9,30,31]. However our estimated effects were larger and 
we speculate these larger rates were due to the severity 
of the Sydney dust storm. Our study showed the dust 
storm period was associated with a decrease in respira- 
tory hospital admissions. While the dust storm was asso- 
ciated with a higher risk of acute respiratory problems, 
these may not have been severe enough to warrant ad- 
mission to hospital. Other studies examining respiratory 
admissions also did not find associations with dust 
storm events [10,28,32]. 

The majority of hospitals in the study area were in- 
cluded. Seasonality was controlled for by smoothing 
functions. We have no data on GP visits. The dust storm 
occurred on a Wednesday, GP clinics would have been 
open during this time. However, we are concerned with 
severe health symptoms and it is likely that GPs would 
have referred patients presenting with severe symptoms 
to an ED. Increased media awareness during the period 
of the dust storm may have biased results, warnings may 
have increased ED visits during this period due to over- 



cautiousness or decreased ED visits due to patients vac- 
ating the study area. We rely on the diagnoses coded in 
the NSW Health datasets, non-differential misclassifica- 
tion is unlikely due to the accuracy of records. 

Public health messages were broadcast in the media in 
response to the dust storm to reduce the incidence of ad- 
verse health outcomes. Coverage included issuing warn- 
ings and recommendations to reduce the likelihood of 
adverse health outcomes: avoid exposure if possible, stay 
indoors (preferably in air conditioned buildings), do not 
exercise, follow asthma plans and seek medical help if re- 
spiratory or cardiovascular symptoms occurred. Because 
the dust storm and consequent public health messages 
had widespread media coverage, the health consequences 
from this dust storm are likely to represent the optimal 
health outcomes that can be hoped for when similar 
events occur in the future. Recent studies evaluating the 
health impact of media alerts on air quality among asth- 
matic patients [33] and the general population [34] sug- 
gest warnings may be associated with behavioural changes 
that reduce exposure to air pollutants. Our study is unable 
to evaluate the efficacy of these warnings; a study involv- 
ing a survey soon after such an event would be the most 
appropriate way to evaluate the impact of media warnings 
[33]. Our study revealed that the dust storm period was 
associated with higher rates in asthma related ED visits. 
We recommend that specific warnings targeting asth- 
matics should be released prior to predicted future dust 
storm events. 

We investigated the effects of the dust storm on the 
broad categories of respiratory and cardiovascular ED 
visits and hospital admissions. Future research could in- 
clude determining the effect of the dust storm on specific 
conditions within these categories such as pneumonia, 
stroke, or myocardial infarctions. Furthermore, utilizing 
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other data sources such as ambulance call-out data may 
complement our results obtained from ED and hospital 
data. 

Conclusions 

The dust storm period was associated with large significant 
increases in asthma ED visits. Respiratory and all-cause ED 
visits also increased, although to a lesser extent. We rec- 
ommend that public health measures including specific 
warnings targeting asthmatics should be implemented dur- 
ing future dust storm events. 
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